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We present the discovery of ferromagnetism in single crystalline Yb(Rho. 73000.27)2812 below Tc ~ 
1.30 K with a spontaneous magnetic moment of 0.1 /^s/Yb along the crystallographic c-axis. This 
is shown by the huge ac-susceptibility peak, the hysteresis in magnetization as well as by the field- 
dependent susceptibility and specific heat data. This discovery motivates the reinvestigation of the 
magnetic order in Yb(Rhi_a;Ooa:)2Si2 with x < 0.27 and in pure YbRh2Si2 under small pressures 
with field along the c-direction to look for low-lying ferromagnetic order. 



PACS numbers: 71.27.+a, 64.70. Tg, 75.50.Cc 

The unambiguous observation of antiferromagnetic 
(AFM) quantum critical points (QCPs) in heavy-fermion 
systems has led to an increasing number of theoretical 
and experimental works in order to understand quan- 
tum phase transitions (QPTs) as deeply as their clas- 
sical counterpart [U |2]- This work is still in progress, 
since, although in few cases the standard spin-density- 
wave model seems to work, in others it does not [31 13] . 
The so-called local QCP scenario has been proposed to 
partially understand some fundamental observations, e.g. 
the oj/T scaling in CeCui_a;Aua; [5], but several other 
features are still unexplained. Nowadays, it is clear that 
many of these quantum critical systems are governed by 
not just one single energy scale and the simple Doniach 
scenario with one single QCP can not be apphed [B]. Pro- 
posals concerning a global phase diagram which incorpo- 
rates the effect of pressure, magnetic field and frustration 
on the magnetic state as well as on the Kondo effect (in- 
cluding Fermi surface reconstruction) have lead to a new 
paradigm [^j^: A system can be tuned across critical 
points of different nature, e.g. with itinerant or local 
character. In this global phase diagram the magnetic 
phase is considered to be antiferromagnetic. However, 
real systems are more complex and often the dynami- 
cal susceptibility is governed by both finite wave vec- 
tor (q=Q) and q=0 critical fluctuations [TUHT^ . Theo- 
retical calculations suggest that when a uniform (q=0) 
clean itinerant ferromagnet is tuned towards its putative 
ferromagnetic (FM) QCP, it becomes inherently unsta- 
ble towards a phase transition of first order or modu- 
lated/textured (q 7^ 0) phases OIH]. Several metals 
have confirmed this prediction [TSHTS] . 

In this respect, the tetragonal heavy-fermion com- 
pound YbRh2Si2 is a prototype system to study, since 
it displays AFM order below a very low temperature 
Tn = 0.07 K [in] and is very close to a peculiar AFM 
QCP [5D] where FM critical fluctuations have been ob- 
served [T0l[2T]. At 2.2 mK another phase transition has 
been observed in magnetization measurements, but its 
origin is still unrevealed [H]. Pressure or Co substitu- 



tion on the Rh site, i.e. Yb(Rhi_a:Coa;)2Si2 , increases 
T/v as well as the strength of the FM fluctuations [23. 24]. 
In addition, a second magnetic transition emerges at 
a temperature Tl < Tn and joins at x = 0.27 
(« 4.3 GPa) [MIES]. Both phases below Tn and Tl were 
believed to be AFM since the application of a magnetic 
field along the magnetic easy axis, i.e. perpendicular to 
the crystallographic c-axis, depresses the transition tem- 
peratures to zero. In this study, we show that the AFM 
ground state of YbRh2Si2 is tuned to a FM one in the 
Co-substituted compound with x = 0.27. This discov- 
ery advises to reexamine the magnetic phase diagram of 
Yb(Rhi_^Co^)2Si2 with x < 0.27 and of YbRh2Si2 un- 
der small pressures to look for a possible FM ground state 
below Tl. 

YbRh2Si2 is strongly anisotropic in terms of the crit- 
ical magnetic field Bn which suppresses the AFM or- 
der at T = 0: Bn,±c = 0.06 T and Bn^\\c = 0.66 T for 
the magnetic field applied perpendicular and parallel to 
the c-axis, respectively 26J. The ordered magnetic mo- 
ment within the AFM phase was determined to be about 
2 • 10~ ^ fiB via muon spin rotation measurements [27] . 
This small value and the very low Tn make it extremely 
difficult to determine the magnetic ordering wave vec- 
tor Q by means of neutron scattering. Therefore, it 
is helpful to increase Tn and the magnetic moment by 
applying hydrostatic pressure, as this usually stabilizes 
magnetism in Yb systems. This has been investigated 
in Refs. [^ [25^50] . Moreover, the second phase transi- 
tion at temperature Tl < Tn emerges above O.IK with 
0.5 GPa. Comparable behavior can be obtained by iso- 
electronic substitution of Co on the Rh site (see Fig. [T^) . 
Good agreement concerning the transition temperatures 
up to about 2 GPa, which corresponds to a Co content 
X — 0.12, was found between the pure system under hy- 
drostatic pressure and the Co-doped samples [211 [13 [31] . 
However, there are some differences: Pressure experi- 
ments have showed that the unit cell c/a-ratio is constant 
up to the highest applied pressures of 21 GPa [29], while 
in Yb(Rhi_2,Co2,)2Si2 an increase of c/a has been seen 
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with decreasing unit cell volume [23] . For small concen- 
trations the change of c/a is linear and amounts to 0.3 % 
for X = 0.27. At a; = 0.27 (marked by the arrow in 
Fig. [T^) , which corresponds to a hydrostatic pressure of 
about 4.3 GPa, and Tl appear to merge to a single 
magnetic transition, while in the pressure experiments 
by Knebel et al. two distinct phase transitions were ob- 
served up to 7 GPa [23]. 

In YbRh2Si2 both FM and AFM fluctuations were 
found to exist [31]. Based on field-dependent specific 
heat and resistivity experiments under pressure, Knebel 
et al. suggested that the AFM state of YbRh2Si2 
changes into a FM one above about 7 GPa [23] . Also in 
Yb(Rho. 73000.27)2812 the increase of chemical pressure 
increases the strength of the FM fluctuations. This has 
been concluded from Curie- Weiss fits performed within 
a fixed temperature range above T/v, where the Weiss 
temperature &w increases with increasing x and even 
switches from negative to positive values at a; = 0.27 [M] . 
However, preliminary investigations suggested an AFM 
state for Yb(Rho. 73000.27)2812, since a magnetic field 
BJ-C of about 0.5 T depresses the transition temperature 
to zero. Application of pressure on Yb systems stabilizes 
its trivalent magnetic state and, therefore, reduces the 
Kondo energy scale. For x = 0.27 a Kondo temperature 
Tk = 7.5 K was found [551 which is about 1/3 of that of 
pure YbRh2Si2. 

In this letter, we present the low temperature 
properties of Yb(Rho. 73000.27)2812 by means of ac- 
susceptibility, specific heat and magnetization in zero 
field and with external magnetic fields applied along the 
crystallographic c-direction. Surprisingly, we find the 
second order phase transition at Tc — 1.30 K to be ferro- 
magnetic. This is shown by the high value of the suscep- 
tibility peak, the hysteresis in magnetization as well as 
by the behavior of the field-dependent susceptibility and 
specific heat. A small spontaneous magnetic moment 
O.l^s/Yb is determined. A second, lower-lying phase 
transition Tj, could not be detected down to the lowest 
achieved temperature of 0.02 K. Such a transition ex- 
ists at lower Go concentration and for pure YbRh2Si2 at 
pressures below 4 GPa (cf. Fig. [I^). Since preliminary 
results suggest that for x < 0.27 the phase below T/y is 
indeed of AFM type, one might find the phase below 
to be FM. 

For all measurements we used the same single crystal 
which was grown by indium flux technique |24j . The size 
of the sample is approximately 1.4mmxl.4mmx0.5mm 
where the last value is the length along the c-direction. 
The cobalt content x = 0.27(1) was measured by energy- 
dispersive x-ray spectroscopy. The low temperature mea- 
surements were performed in Oxford Instruments '^He- 
He-dilution refrigerators in a temperature range 0.02 < 
T < 5K and magnetic fields up to IT. The magneti- 
zation was measured with a Faraday magnetometer |33| 
and the ac-susceptibility with a standard susceptometer 
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Figure 1: (Color online) (a) Excerpt of the magnetic phase 
diagram of Yb(Rho. 73000.27)2812 from Ref. [21]. Open and 
closed circles correspond to the transition temperatures Tl 
and Tjv, respectively. For x < 0.12, the data were taken 
from susceptibility measurements 25 , while for larger Oo 
content the data points were obtained by specific heat mea- 
surements [21]. The red arrow marks the position of the 
present system, where we observe FM order. The lower 
axis marks the corresponding hydrostatic pressure pj; on 
YbRh2Si2 (cf. Ref. f24l). (b) Ac-susceptibility x {T) of 
Yb(Rho. 73000.27)2812 in different magnetic fields with B\\c. 
At zero field a sharp peak is observed at Tc = 1.30 K. With 
increasing magnetic field the peak decreases rapidly and its 
position shifts towards higher temperatures. Inset: Temper- 
ature dependence of the imaginary part of the susceptibility 
X"{T). (c) Ac-susceptibility x'{T) of Yb(Rho.73Coo.27)2Si2 
parallel (black) and perpendicular (gray) to the c-direction in 
zero magnetic field. A magnetic field _B_Lc = 0.5 T (violet) 
suppresses the phase transition towards lower temperatures 
(marked by arrow). 



with a modulation field Brms = 11 /^T. The absolute val- 
ues were obtained by matching the data to measurements 
done in a commercial Quantum Design SQUID Vibrating 
Sample Magnetometer. The specific heat was measured 
with a semi-adiabatic heat pulse technique [Mj and a 
Quantum Design Physical Properties Measurement Sys- 
tem. 

Figure shows the ac-susceptibility as a function of 
temperature in various magnetic fields. At zero field 
a sharp maximum is observed at Tc = 1.30 K with a 
huge absolute value of 73 • 10~^ m'^/mol. With the molar 
volume V^oi = 46.7 • lO'^m^/mol, we obtain the di- 
mensionless volume susceptibility Xvoi — x' /^moi — 1-5 
which corresponds to a lower limit of the intrinsic sus- 
ceptibility because of demagnetization effects due to the 
sample shape. However, the large value of Xvoi points 
to a FM phase transition, or at least a phase transition 
with a FM component along the c-axis. The applica- 
tion of an external static magnetic field strongly influ- 
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Figure 2: (Color online) Specific heat C/T vs T of 
Yb(Rho. 73000.27)2812 in different magnetic fields. At B = 0, 
C/T features a mean- field-like anomaly at Tc ~ 1.30 K. In- 
creasing B, the anomaly is shifted towards higher tempera- 
tures. Inset: The black curve shows the specific heat in zero 
field down to 0.03 K. The red curve shows the same data with 
the nuclear Schottky contribution subtracted. 



Figure 3: (Color online) Isothermal magnetization curves for 
different temperatures with B\\c. The field axis is shifted by 
0.008 T in order to compensate for the remanent field due to 
the magnet gradient field. The arrows indicate the field sweep 
direction. The spontaneous magnetization is very small, and 
the obtained total signal is close the magnetometer's resolu- 
tion limit. 



ences the size of the peak. A tiny field B ~ 0.001 T 
reduces its size by more than a factor of 2. Increasing 
B, the maximum shifts towards higher temperatures as 
expected for the crossover temperature (which separates 
field-polarized and paramagnetic state) of a ferromagnet. 
The inset of figure [T|3 shows the imaginary part x"(T) of 
the susceptibility. x"{T) exhibits a clear feature within 
the critical region between 1.15 and 1.35 K at very low 
fields indicating the presence of dissipation (non-zero hys- 
teresis), characteristic of ferromagnets. At fields above 
0.01 T the imaginary part vanishes. The susceptibility of 
Yb(Rho. 73000.27)2812 is reminiscent to the situation in 
YbNi4P2 which is located close to a FM-QCP [35]. In- 
terestingly, in both systems the magnetic hard direction, 
which has a tiny susceptibility at enhanced temperatures, 
eventually becomes much larger slightly above Tc- This 
can be seen by comparing our data with the suscepti- 
bility for BJ-C, which is shown in figure [l]:. Moreover, 
x'{T) for BJ-c — 0.5 T demonstrates the suppression of 
the phase transition towards lower temperatures unlike 
the crossover for the B\\c case. The phase diagram in 
terms of specific heat, magnetization, ac-susceptibility 
and electrical resistivity for BJ-c will be published in a 
forthcoming paper [36 . 

The specific heat of Yb(Rho. 73000.27)2812 is depicted 
in Fig. [I At zero field, C{T)/T shows a mean-field-like 
transition at Tc = 1.30 K in agreement with the peak 
temperature found in }('{T). The application of a mag- 
netic field shifts Tc towards higher temperatures and the 
maximum broadens as expected for a FM phase transi- 



tion. In the inset of Fig. [2] the low temperature specific 
heat (black curve) is shown at i? = 0. The subtraction of 
the quadrupolar contribution of the Yb nucleus gives only 
a tiny correction to the data [29l [37] . A more important 
contribution can be ascribed to the magnetic moment of 
the 4/-electrons which creates a magnetic field at the Yb 
nucleus. This leads to an additional Schottky splitting of 
the nuclear dipole moment which gives rise to a con- 
tribution to C{T)/T even at zero field. We suppose that 
the electronic contribution becomes constant at lowest 
temperatures and that the low temperature divergence 
is due to nuclear Schottky splitting induced by the spon- 
taneous magnetization. In order to obtain a constant 
C{T)/T, we subtracted the nuclear Schottky contribu- 
tion from our specific heat data by adjusting the spon- 
taneous moment of the 4/-electrons. For this we used 
the hyperfine coupling constant determined from NMR 
measurements [38j and the size of the nuclear dipole mo- 
ment [37]. The result is depicted by the red curve in the 
inset of Fig. [2] which exhibits a significant deviation from 
the raw data below 0.2 K. From this procedure we can 
extract a 4/ magnetic moment of 0.20(5) /is/Yb. 

Figure [3] illustrates the magnetization M{B) as a func- 
tion of magnetic field at various temperatures. At the 
lowest achieved temperature of 0.08 K a clear hystere- 
sis is observed with a coercive field of about 0.08 T and 
a tiny spontaneous magnetic moment 0.05/iB/Yb. The 
hysteresis shrinks with increasing temperature and dis- 
appears for T > Tc- The size of the magnetic moment 
is smaller than the value from the specific heat measure- 
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Figure 4: (Color online) Magnetization as a function of mag- 
netic field B\\c at T = 1.8K. The Van Vleck contribution 
Xvv = 0.172 • 10~^m^/mol was subtracted from the raw 
data, leaving a saturation magnetization at 7T of about 
0.11/^fl/Yb. Inset: data of Xdc ~ M/B between 2 and 
300 K. Below 20 K a Curie- Weiss fit (red curve) was carried 
out as described in the main text. 



mcnts. This might be due to a complex magnetic struc- 
ture (e.g., spiral) with the total 4/ magnetic moment 
being larger than its component in c-direction. 

In Fig. |4] the magnetization M is plotted as a function 
of B above Tc at T = 1.8 K (black curve) . Up to a field of 
7 T the curve does not saturate which is due to a consid- 
erable Van Vleck contribution. To determine this contri- 
bution we performed a measurement of xdc = M/ B as a 
function of T in an external magnetic field of 0.5 T. The 
data are shown in the inset of Fig. |4] The broad shoul- 
der above 20 K corresponds to the first excited doublet 
of the J = 7/2 multiplet of the Yb'^+ atom split by the 
crystalline electric field. We have fitted the region below 
20 K with the function x = C/{T — Qw) + Xvv where 
C is Curie's constant, Qw = 1-24 K is the Weiss tem- 
perature and Xvv = 0.172 • 10~^ m^/mol is the constant 
Van Vleck term. From the Curie's constant we can calcu- 
late the corresponding magnetic moment 0.22 hb along 
the c-direction. This value is in nice agreement with the 
saturation magnetic moment of 0.25(5) /x^/Yb which is 
expected for the doublet ground state from electron-spin- 
resonance (ESR) results on Yb(Rho.73Coo. 27)2812- There, 
a (7-factor oi g = 0.5(1) was determined at T = 5 K along 
the c-direction [33] . In the main panel of Fig. |4] we have 
subtracted the Van Vleck contribution xvv from M{B) 
and obtained the blue curve which seems to saturate 
above 7T at about O.ll/is/Yb. This value is substan- 
tially smaller than the values obtained from the Curie- 
Weiss fit and the ESR results, respectively. We ascribe 
this to the Kondo screening which still exists at T = 1.8 K 



(< Tk = 7.5K) and 7T (< Bk = {kB/gfiB)TK ~ 22T). 

To summarize, we have shown that 
Yb(Rho.73Coo. 27)2812 is a ferromagnet with a Curie 
temperature Tq — 1.30 K and a remanent magnetization 
along the crystallographic c-direction. This is evidenced 
by the huge absolute values of the ac-susceptibility, the 
hysteresis in the magnetization isotherms at T < Tc, 
and by the evolution of the transition in magnetic field. 
A spontaneous magnetic moment of 0.05/XB/Yb has 
been measured along the c-direction in the FM phase, 
which is smaller than the value estimated from specific 
heat (0.20(5) /Lts/Yb) which might be due to a complex 
magnetic structure. Neutron diffraction experiments 
were also performed on Yb(Rho.73Coo. 27)2812, but, so 
far, no magnetic structure could be resolved below Tc- 
As a matter of fact, this is expected for a magnetic 
ordered moment of the order of 0.2 /xb, which is below 
the resolution limit of common neutron diffractometers. 
The negative neutron diffraction results also suggest 
that, if an AFM component is present at all, it must 
be in the same order of magnitude. Above Tq, we 
recognized the c-direction to be the magnetic hard axis 
and determined a magnetic moment of 0.22 /is from a 
Curie- Weiss fit below 20 K. After subtracting the Van 
Vleck contribution, the magnetization results taken at 
1.8K reveal O.ll/is/Yb at 7T, which proves that the 
Kondo effect is operating. Our discovery opens the 
question whether the low-lying phase below which is 
found for small Co concentrations in Yb(Rhi_a:Coa;)28i2 
or in the pure YbRh2Si2 (under pressure) is ferromag- 
netic, too. This might imply a field-induced FM QCP 
at Bl < Bn- 

We acknowledge fruitful discussions with O. 8tockert, 
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